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ABSTRACT 



The Maximum Likelihood estimators for the Neyman Type A 
distribution parameters are very difficult to compute. In 
this thesis, the Empirical Probability Generating Function 
is used to provide estimators that are easier to compute and 
have asymptotic efficiency at least as high as 97% of that 
for the Maximum Likelihood estimators over most of the param- 
eter space considered. The estimators found by this method 
are consistently better than the Method of Moments and the 
Method of Zero Frequency estimators with respect to asymptotic 
efficiency. The considerations of preference in using one 
method over another are discussed. 
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I. 



INTRODUCTION 



Several 'contagious' distributions, including the Neyman 
Type A distribution, were derived by Neyman [1] in the course 
of modeling the reproduction behavior of certain types of 
plants and bacteria. More currently, these models come under 
the heading of branching processes. Several methods of pa- 
rameter estimation have been proposed for the Neyman Type A 
distribution with varying degrees of success. The two most 
frequently used are the Method of Moments and the Method of 
Zero Frequency. The performance of an estimation scheme can 
be measured by the asymptotic efficiency of its estimators 
as defined by Wilks [2] . This paper proposes a method for 
estimating the parameters of the Neyman Type A distribution 
by using the Empirical Probability Generating Function. 

For the Neyman Type A distribution, the Method of Moments 
provides estimates that are quick and easy to compute, al- 
though the efficiency of these estimators is not always high. 
It is well known that Maximum Likelihood estimators are 
asymptotically efficient, but they are very difficult to com- 
pute for the Neyman Type A distribution. 

The method of estimation by using the Empirical Probabil- 
ity Generating Function provides estimators with asymptotic 
efficiency close to that of the Maximum Likelihood estimators 
over most of the parameter space considered. Extensive tables 
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are presented along with an algorithm that is relatively easy 
to use on current programmable calculators. 

The estimators found by applying this method have con- 
sistently higher asymptotic efficiency than the Method of 
Moments and the Method of Zero Frequency estimators. However, 
larger sample sizes may be required to cause the algorithm to 
converge on a unique pair of estimates . 

A brief outline of the report is presented as follows. A 
derivation of the Neyman Type A distribution with a list of 
some of its mathematical properties is contained in section 
II. The equations for the Maximum Likelihood estimators and 
asymptotic efficiency appear in section III. The mathemat- 
ical details of using the Empirical Probability Generating 
Function to provide estimators are found in section IV, along 
with a pictorial display of asymptotic efficiency for these 
estimators. An algorithm for computing the estimates of this 
method is presented in section V. The method of using the 
Empirical Probability Generating Function is compared to the 
Method of Moments in Section VI. The algorithm described in 
section V is illustrated by a numerical example along with 
computer studies in section VII. Tables for use with the 
above algorithm are found in section VIII and a summary of 
the results appears in Section IX. 
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II. NEYMAN TYPE A DISTRIBUTION 



Neyman's Type A distribution can be used to model decay 
and reproduction processes which fit the assumptions of com- 
pound Poisson processes. A derivation of this distribution 
provides an understanding of its structure, usefulness, and 
properties. For these reasons it is presented in the fol- 
lowing form. 

Consider a decay process in which radioactive nuclei 
decay to produce daughter nuclei according to a Poisson proc- 
ess with rate X. Also suppose that the number of daughter 
nuclei produced after each decay are independent random 
variables having a common Poisson distribution with parameter 
rr^. Let N(t) = the number of decays in the time interval 

(0 , t] . 

Y^ = the number of daughters produced from the ith 
decay . 

X = the total number of first generation daughters 
produced in the time interval ( 0 , t ] . 

Then 



N(t) 

X = Z Y. . 
i=l 

Because N(t) is a Poisson process with rate X, 

-Xt. (Xt) n 

P [N (t) = n] = e nl for n = 0,1, 
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Conditioning on N(t) is used now to find the probability mass 
function of the random variable X. For fixed N(t), X is a 
Poisson random variable with parameter nir^ when N(t) = n. 



Thus 



-nrr^ (nn^) 



x 



P (X = x I N(t) = n) = e x! for x = 1,2,3, 

The event that no daughters are produced in the first gener- 
ation may be described by two cases, namely 

1) X = 0 if N(t) = 0 or, 

2) X = 0 if N ( t) 1 0. 

From case 1) P(X = 0 | N(t) =0) = 1 is obvious. 

-nm_ 

And from case 2) P (X = 0 | N(t) = n f 0) = e 

Now, using the total probability theorem. 



P = P(X = 0) = Z P(X = 0 | N (t) = n) PCN(t) = n) 

° n=0 



or 



-Xt 



-run. 



-Xt (Xt) 



n 



= e 



Z 

n=l 



n! 



-m. 



-Xt(l - e 



) 



= e 



It is also obvious that for X=x>0, P(X=x | N(t) = 0) 
And by using the same conditioning argument as above it is 
evident that. 



= 0 . 
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/ \ X 

rim,-* ( nnu ) > , / . \ ri 

<» 2 2 -At (At) 

P = P (X = x) = Z e x! e nl 

x n=l 

These two formulas may be combined into one expression which 

is given in Shenton [3] as 

-m, „ _ 

P.. = e 1 m n x [ 0 X + S1 X + 6 2 X + 6 3 X + ] (2.1) 

X £ T » I ^ I 



where 



-nu 



m^ = At 



6 = m^e 



and 



1 if x = 0 
0 if x ? 0 



This probability distribution was originally derived by 
Neyman and is called the Neyman Type A distribution. Other 
forms of equation (2.1) may be found in Johnson and Kotz [4]. 

Some of the properties of this distribution are needed in 
the further development of this paper and are 

1. Probability generating function, 

xtu (t-1) 

X -m, [ 1 - e ] 

E ( t ) = e 1 (2.2) 

2. Moment generating function, 



uX 



E (e ) = e 

By using the well-known formula 

uX 



-m^ [ 1 - e 



m 2 (e 



u 



- 1 ) 



(2.3) 



a" E(e ) 



du 



u = 0 



= E (X / 



the following are obtained 
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- 




3. 



E(X) 



(2.4) 



m l m 2 



4. E(X 2 ) 

5. E(X 3 ) 

6. VAR ( X ) 



m^ir^ ( 1 + + m^rr^ ) 

2 

m l m 2 ( 1 + 3m 2 + 3m l m 2 + 3m i rn 2 + 

2 2 2 , 

m 2 + m l m 2 

(1 + m 2 ) 



(2.5) 



( 2 . 6 ) 

(2.7) 
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III. EFFICIENCY AND MAXIMUM LIKELIHOOD EQUATIONS 



A method for estimating the parameters of a distribution 
should yield estimators with high asymptotic efficiency and 
it is well known that Maximum Likelihood estimators have 
asymptotic efficiency equal to 1. However, the Maximum Like- 
lihood equations may be difficult to solve and such is the 
case for the Neyman Type A distribution. Replacing some of 
these equations in the Maximum Likelihood system may produce 
a system of estimating equations that is more readily solved 
and still produces estimators with high asymptotic efficiency. 
Estimating equations considered here are those that equate a 
statistic with its expected value. This paper presents one 
such system for the Neyman Type A distribution. 

It can be shown (see, for example. Read [5]) that the 
asymptotic efficiency of the estimators of such a substitute 
system is 

EFF = ^ 

I A | 

where | A | = information determinant of the probability dis- 
tribution considered and 

/\ yX T 

M = lim [ nE(0 - 0) (0 - 0) ] 

n-><» 

with 0 = vector of distribution parameters 
and 0 = vector of estimators of 0. 
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,e expression for EFF can be simplified if 
vector of estimating equations 



A = the matrix E f|^ j and C = lim [ nE(hh T ) ] 

' n-*- 00 




EFF = . (3.1) 

I A | | C | 

5 Maximum Likelihood equations for the Neyman Type A 
htion are derived in [3] and the following forms are 
i [4] 

A A 

x = m^rr^ 
n P . 

2 (x. + 1) x i = nx (3.2) 

i=l 1 P 

x . 

l 

= the sample mean 
= the sample size. 

: should be noted that equation 2) is an extremely corn- 
id function of m^ and m 2 • To compute the efficiency 
:,j other set of equations as substitutes for the equa- 
.n (3.2), it is necessary to compute the information 
Lnant of the Neyman Type A distribution. The informa- 
sterminant has been derived in [3] and there it is 
that 

2 

= 1 [ (1 + m 2 )Q ~ m ]_ m 2 + m i m 2 + m 2^ ^ (3.3) 

m l m 2 



1 ) 



Le 



) 
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IV . ESTIMATING EQUATIONS 



The Maximum Likelihood equations are quite difficult to 
solve for the Neyman Type A distribution as seen from the 
equations in (3.2). Namely, the second equation is the one 
that produces the complication. Replacing the second equation 
in (3.2) with a different one will allow for simpler calcu- 
lations if this new equation is chosen judiciously. However, 
the idea of maintaining high efficiency should also be con- 
sidered when making this choice of a substitute. The use of 
the Empirical Probability Generating Function (hereafter 
referred to as EPGF) seems to be appropriate in attaining 
this goal. This was suggested by the successful application 
of this choice in the case of the Negative Binomial distri- 
bution as demonstrated by Caglayan [7] . Since the EPGF 
contains the independent variable t, the solution of the new 
system of estimating equations can be "tuned" so as to achieve 
the estimators of maximum efficiency capable from the new 
system. Considering the form of the probability generating 
function for the Neyman Type A distribution, it is easily 
amendable to programming on current programmable calculators. 
With this motivation in mind, the first Maximum Likelihood 
equation in (3.2) and the EPGF equation which follows are 
considered as a possible alternative to estimating the pa- 
rameters for the Neyman Type A distribution. 
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The EPGF is defined as 



n x . 

EPGF = 1 E t 1 

n i=l 

where n = the sample size and = the ith sample value. The 
expected value of the EPGF is the probability generating func- 
tion of the distribution of the random variable X and is given 
by equation (2.2). The system of estimating equations con- 
sidered above is then 

1) x = m-^m 2 (4.1) 

rt»2 ( t-1) 

n x. -m. [ 1 - e ] 

2) 1 E t 1 = e 1 (4.2) 

n i=l 

The notation m^* and m 2 * is used to denote the estimators 
of m^ and m 2 found by solving equations (4.1) and (4.2). To 
solve these equations the variable t must have a value in the 
interval [0,1) to insure that the sign of the left-hand side 
of equation (4.2) is compatible with the right-hand side. 

The value chosen for t (which will be denoted as t*) is found 
by maximizing the asymptotic efficiency of the estimators m^* 
and m 2 *. 

To compute the asymptotic efficiency of m^* and m 2 * 
equation (3.1) is used. To reduce the amount of typing in 
the following expressions let 

nu (t-1) 

-m, [ 1 - e ] 

GCt) = e 1 
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The vector of equations (h) is given by 



x m-j in^ 



= 1 



The matrix A is 



n 


X. 

1 


Z 

i=l 


t 




31^ 


3m^ 


3m 2 


3h 2 


3h 2 


3m^ 


3m 2 



from which it follows that its elements are 

A, 



11 


— 


m 2 




(4.3) 


12 


= 


" m l 




(4.4) 


21 


= 


m 2 (t-1) 

1 - e 


G ( t ) 


(4.5) 


22 


m„ (t-1) 

= m 1 (l-t)e G ( t ) . 


(4.6) 



The matrix C is 



lim nE 
n->°° 



h 2 h l 



h l h 2 



and it follows that its elements are 



= lim nVAR(x) = m^i^d + m 2 ) 
n-*-°° 



(4.7) 



n x . x . 

C,. = lim nCOV (x, 1 Z t X ) = COV( x., t 1 ) 

n-*-°° n i=l 



'12 



the last step resulting from x^ and x ^ being independent when 

i 7* j • 
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Continuing, 



x. x . 

C 12 = E(x i t } " E(,x i )E(t ± ) 

and using the fact that for well-behaved distributions 
eluding Neyman's Type A) 

X . x . 

E(x.t 1 ) = tE [ _d(t X ) ] = tG’ (t) 

dt 

then 



m 2 (t-1) 

C 12 = m l m 2 G ^ [ te - 1 ]. 



By symmetry = C 21* 



n x . x . 

= lim nVAR(l E t 1 ) = VAR ( t 1 ) 

n-*- 00 n i=l 



thus 



C 22 = G(fc2) " [GCt)]2 - 
Using equations (4.3) - (4.6) 

rtu (t-1) 

| A| = m 1 G(t) [ 1 + (m 2 (t-l) - l)e z ] 

and using equations (4.7) - (4.10) with 

m 2 ( t— 1 ) 2 

B = 1 + m 2 + m-^m 0 (te - 1) 

then 

| C j = m 1 m 2 [ (l+m 2 )G(t 2 ) - G 2 (t) B ] . 

Now substituting (4.11) and (4.12) into (3.1) results 

2 m 2 ( t— 1 ) 2 

m.G z (t) [ 1 + (m_t - m„ - 1) e ] 

eff = — = 

|A|m 2 { (l+m 2 )G(t 2 ) - G 2 (t) B ] 



(in- 

(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 
in 

(4.13) 
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t* is the value of t used to solve equations (4.1) and (4.2). 
It is defined as 

EFF ( t* ) = max EFF(t). 

0<t<l 

In order to use this definition the values of m^ and m 2 

must be given. t* and EFF(t*) have been computed for known 

m^ and m 2 and are presented in Table 1. An algorithm for 
using these tabled values to find the estimates m^* and m 2 * 
is presented in section V.B. 

To describe the relationship between t* , m^ and m 2 , a 
three-dimensional plot was made (Fig. 1) for m^ and m 2 in the 
range .01 (.01). 09, .1(.1).9, 1(1)9. Fig. 2 is a corresponding 
plot of EFF(t*) for the same set of values for m^ and • 

Figures 1 and 2 are representations of the values found in 

Table 1. 

For given values of m^ and m 0 , t* is the value of t which 
maximizes the efficiency function EFF(t). The values of t* 
are not necessarily continuous as m^ and m 2 are varied. The 
efficiency function determined by equation (4.13) is bimodal 
for certain values of m^ and m 2 and these modes occur near 
the endpoints of the interval [0,1). The global maximum of 
EFF ( t) on this interval may occur at either mode for different 
values of m^ and m 2 , which is the cause for the discontinui- 
ties of t* as seen in Fig. 1. The jump in the surface occurs 
when m^ is about 3.0 for all values of m 2 > 2.0. Further, it 
is interesting to note that t* is approximately zero for 
m^ < 1.0 and all values of m 2 , except for a singular point at 
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m^ = .09 and m 2 = .01. EFF(.t*) appears to be a continuous 
function of and m 2 as shown in Fig. 2 and is very close 
to 1.0 over much of the range for m^ and m 2 described above. 

The method used for determining t* is discussed in 
section V.A. 
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1.0 




Fig. 1 t* vs . parameters and m 2 
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Fig. 2 EFF(t*) vs. parameters m^ and m 2 
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V. CALCULATING TECHNIQUES 



Estimation of the parameters m^ and m 2 of the Neyman Type 
A distribution using the method outlined in section IV in- 
volves solving a nonlinear programming problem; specifically, 
finding a solution (m^*, ir^*, t*) that maximizes equation 
(4.13) subject to the constraint equations (4.1) and (4.2). 
Since this requires a large amount of computing power, an 
approximating sequential method is proposed that can be per- 
formed on a calculator. In order to do this the values of t* 
must be known for any m^ and m 2 pair. These t* values are 
presented in Table 1 along with the associated maximum EFF 
for a wide range of values of m^ and m 2 . Bivariate linear 
interpolation as discussed in Abramowitz and Stegun [8] can 
be used to approximate t* between the values of m 1 and m 2 
given in Table 1. 

A. METHOD OF COMPUTING t* FOR TABLE 1 

Table 1 was prepared by fixing m^ and m 2 at constant 
values in the range .01 to 9.0 and then conducting a single 
variable search on the function EFF (as a function of t only) . 
A Golden section search procedure was employed that terminated 
at the value of t* equal to the midpoint of the searched in- 
terval on which the value of EFF(t) was increasing but did not 
change by more than 10 across the interval. Because of the 
bimodality of EFF(t) on some regions, a search was initiated 
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from the right end of the interval [0,1) and a separate search 
started from the left end. The global maximum was then found 
by comparing the resulting local maximums . 

B. METHOD OF SOLVING FOR n^* AND m 2 * 

Given a sample of n observations on the random variable X 
assumed to have an underlying Neyman Type A distribution of 
which x^ represents the ith sample observation, the algorithm 
for solving equations (4.1) and (4.2) to find the estimators 
m^* and m 2 * of the parameters m^ and m 2 consists of the fol- 
lowing steps : 

STEP 1 - Compute the Method of Moments estimators from 

the following equations 

n 

x = 1 I x . 

n i=l 1 

? n — 2 

S^ = 1 Z (x. - x) z (5.1) 

n-1 i=l 1 

m^ = x/m 2 

m 2 = - 1 

X 

STEP 2 - Find t* in Table 1 corresponding to m^ and m 2 . 

Note: If t* = 0. and the sample contains no zero values, set 

m^* = m^ and m 2 * = m 2 then STOP. 

STEP 3 - Solve the system of equations (4.1) and (4.2) 

for m^* and m 2 * with t = t* by using a Newton-Raphson pro- 
cedure. Solving equation (4.1) for m^ and substituting this 
into equation (4.2) yields 
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x(g-l) - In ( t X ) 



m. 



where 



(5.2) 



m 2 (t-1) 

g = e 

Applying the Newton-Raphson method on equation (5.2) produces 
as the (k+l)st iterate of m 2 * 



1 - e 



m k (t-1) 



m. 



Ck+1) 



* m k t 1 + 



+ m k In (t ) 



1 + e 



m k (t-1) 



( m k t - - 1) 



where for ease of typing m, = m 9 and iru = m 0 (the Method 

K z (k) z (l) z 



of Moments estimator) . When termination of the above iteration 

scheme occurs by the appropriate criteria as set by the user, 

set m~* = iru and compute m. * from equation (5.1) . 

(n) 1 

STEP 4 - If |m 2 * - m 2 | < e STOP, where £ is set by the 

user. Otherwise return to STEP 2 after replacing fn 2 with m 2 *. 

The convergence of this algorithm depends on the sample 
size n and how well the Method of Moments estimators approxi- 
mate the true values of m^ and m 2 . Experience has shown that 
cycling between estimates can occur when the m^* and m 2 * 
values are close to the points where t* is discontinuous as 
determined from Fig. 1. This behavior was observed when the 
above algorithm was applied to data of sample sizes as large 
as 50. For most cases the estimates obtained by applying the 
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algorithm converge after two or three iterations. It should 
be noted from STEP 2 that if t* = 0 and the sample contains 
no zero values, then the Method of Moments estimators are the 
default values obtained when using the algorithm. 

C. METHOD OF DETERMINING Q 

The value of Q as defined in equation (3.4) was computed 
according to equation (3.5) for the purpose of determining the 
asymptotic efficiency of the EPGF estimators given in Table 1. 
The xth power moment of a Poisson distributed random variable 
(denoted by s in equation (3.5)) was computed by finding the 
maximum term of the series and using this as a scaling factor 
to keep the partial sums within the limit of an IBM 360/67 
computer, which was used for all computations. Actual compu- 
tations were performed with 16 digit precision using logarithms. 

Truncation of the series occurred when the first scaled term 

""161 1 

after the maximum term was found to be less than e ’ . The 

value of Q was calculated by computing the successive partial 

sums of q until the relative error between the last two was 
^x 

_ 7 

less than 10 (the last term was not added to the partial sum 
computed previous to it), i.e., when 

q m+l 7 

< 10 . 

m 

Z q x 
x=0 

A discussion of this stopping rule for the particular appli- 
cation described above may be found in Katti and Gurland [9] . 

As seen in Table 1 there are a few instances where the value 
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of EFF(t*) was computed as a value greater than 1.0, which is 
impossible. It is felt by this author that in these occur- 
rences the value of Q was not computed accurately enough. 

The reason for this is that when m^ or m 2 was large over 300 

terms were required in the partial sum of Q in order to reach 

the stopping rule criterion given above. However, since the 
size and frequency of these inconsistencies were small, it is 
felt that the other results presented in Table 1 are not com- 
promised. This is justified by noting the continuous behavior 
of EFF(t) as given by equation (4.13) and observing the few 
values of EFF(t*) > 1.0 in Fig. 2. It is not known why these 

occurrences did not happen at the extreme values of m-^ = 9.0 

and m 2 = .01. As shown in Fig. 2, they appear to be isolated 
points of irregularity in the computation of EFF(t). 
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VI. COMPARISON WITH ESTIMATION BY THE METHOD OF MOMENTS 



To determine the utility of the EPGF method of estimation 
for the Neyman Type A distribution, a comparison with Method 
of Moments estimation was made. This comparison was performed 
by computing the ratio of the asymptotic efficiencies of the 
two schemes. The estimators of m 1 and m 2 by the Method of 
Moments are found by solving the equations of (5.2) and is 
actually done in STEP 1 of the sequential algorithm presented 
in section V.B. The EPGF method is an attempt to improve on 
the Method of Moments estimation of m^ and m 2 . The amount of 
improvement is determined by the increase in asymptotic ef- 
ficiency gained by applying the EPGF method relative to that 
of the Method of Moments. In order to find this increase, 
the asymptotic efficiency of the Method of Moments must be 
found for the Neyman Type A distribution. Using equation 
(3.1), this was done and is shown below. 

The determinants of A and C for the Method of Moments are 
determined in a manner similar to section IV. By using the 
equations in (5.2) the following results can be obtained: 

m ^ 

= -m x 

2 

= • nr 2 

= “ m ]_ ~ 2 m -^ m 2 



i n 



12 



21 



22 
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thus 

2 

| A | = m 1 m 2 . (6.1) 

Also 

C 1L = VAR (X ± ) = rn^Cl + m 2 ) 

C 12 = E(X i 3) " 3E(X i 2 )E(X j ) + 2 [ E(X i ) ] 3 

2 

C 12 = + 3m 2 + m 2 ) 

and 

C 21 = C 12 * 

From formulas derived in [5] , 



C 22 = E(X i - X ) 4 - [ VAR(X i ) ] 2 

2 2 3 

C 22 = m l m 2 ^ 2m^m 2 (1 + m 2 ) + (1 + 7m 2 + 6m 2 + m 2 ) ] 

so that 

| C | = m^ m 2 ( 2m^m 2 (1 + m 2 ) + 2m 2 + 2m 2 + m 2 1* (6*2) 

Substituting equations (6.1) and (6.2) into equation (3.1) 
yields 



m. 



EFF 



mm 



A | [ 2m^m 2 (1 + m 2 ) 3 + 2m 2 + 2m 2 2 + m 2 3 ] 



(6.3) 



where EFF^ represents asymptotic efficiency of the Method of 

Moments estimators m, and iiL . The ratio EFF/EFF formed from 

12 'mm 

equations (4.13) and (6.3) represents the relative merit of 
the EPGF to Method of Moments estimation. 
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Values of t* and EFF/EFF are tabulated in Table 2 for 

mm 

the same range as covered by Table 1. This comparison of 
the two methods is most readily observed in Fig. 3, which is 
a three-dimensional representation of Table 2. The EPGF 
method is better than the Method of Moments in producing 
estimators of high asymptotic efficiency for the range of 
parameter values covered in Table 2. This is clearly demon- 
strated by Fig. 3. The EPGF method reduces to the method of 
estimation by zero frequency when t* = 0. In this case the 
contours of Fig. 3 can be found in [4] and [9]. 
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Fig. 3 Efficiency of the EPGF method relative 
to the Method of Moments vs . the 
parameters m^ and m 2 
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VII. NUMERICAL EXAMPLES 



EXAMPLE 1 

To illustrate the EPGF method and in particular the al- 
gorithm of section V.B, the data of insurance claims caused 
by accident or sickness cited by Shenton [3] are used and 
reproduced below. The data are assumed to be observations 
from a Neyman Type A distribution and three methods of esti- 
mating the parameters m^ and m2 are considered; 

1) the Method of Moments, 

2) the EPGF method, and 

3) the Method of Maximum Likelihood. 

The data are : 



Value 



Frequency 



0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 



187 



185 

200 

164 

107 



68 

49 

39 

21 

12 

11 

2 

5 

2 

3 

1 



sample size 



1056 



sample mean 



2.8059 



sample variance = 6.4106 
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When the EPGF algorithm of section V.B is applied to this 
data, the following results are obtained: 

a) from STEP 1 the Method of Moments estimators are 
computed to be 

m^ = 2.184 and ir^ = 1.285 

b) from STEP 2 the t* value corresponding to m^ and rr^ 
is found by bivariate linear interpolation in Table 1 as 
follows ; 

for m^ = 2.0 and m 2 = 1.285 linear interpolation gives 
t* = .0748, 

for m^ = 3.0 and m 2 = 1.285 linear interpolation gives 
t* = .3774, and 

for m^ = 2.184 and m 2 =1.285 linear interpolation gives 
t* = .1305 

c) from STEP 3 the solution of (5.2) yields 

m^* = 2.638 and m 2 * = 1.064 

where termination of the Newton-Raphson method was set for an 
absolute error of .0001 

d) from STEP 4, where c = .01, a second application of 
the algorithm is required and the above m^* and m 2 * is used 
as the new starting point in STEP 1. This second iteration 
then yields 

t* = .2736, m L * = 2.609, and m 2 * = 1.076. 

Since the termination criteria of STEP 4 is not met a third 
iteration is performed resulting in the following 
t* = .2653, n^* = 2.611, and m 2 * = 1.074 
with the termination criteria being satisfied. 
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The computation of Maximum Likelihood estimators for the 
Neyman Type A distribution is very complicated and only ap- 
proximate Maximum Likelihood estimates for the above data 
are given in [3]. To get an idea of the computational effort 
involved, these estimators are only the first iteration re- 
sults and required several hours on a desk calculator in 
] 949 . 

Maximum Likelihood estimators for the above data were 
found using a TI-59 programmable calculator and the method 
presented in [6]. The range of observed values was small 
enough so that data storage was not exceeded on the calculator. 
The following results were obtained in approximately thirty 



minutes of computing time, 

/\ /v 

Iteration m l m 2 



1 2.475 1.134 

2 2.515 1.116 

3 2.516 1.115 

4 2.516 1.115. 



In comparison, the computing time used to perform the 
EPGF method was approximately five minutes and the data 
storage requirement consisted of only keeping eight numbers, 
which were not influenced by the data values. 

EXAMPLE 2 

Computer simulation was used to generate eighty-one data 
sets of random numbers from the Neyman Type A distribution. 
Computer program 4 in Appendix A was written to simulate the 
model presented in section II and the data sets produced are 
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given in Appendix B. This program uses a random number gen- 
erator that provides random deviates from a Poisson distri- 
bution and is a local subroutine. For each data set produced, 
the sample size was fifty and the range of values used for m^ 
and m 2 was 1(1)9. 

The EPGF estimates for each data set were found by apply- 
ing the algorithm presented in section V.B. The t* value 
required in STEP 2 was computed by actually maximizing the 
efficiency equation (4.13) with the estimated values of m^ 
and n^. The results of these computations are found in Table 
3. The termination criteria of STEP 4 was set at £ = .0001, 
so that more iterations were required than would normally be 
expected. 

In only two of the eighty-one cases considered, no im- 
provement was achieved over the Method of Moments estimators. 
These two occurrences are noted in Table 3. For the case 
m^ = 4.0 and m 2 = 7.0, the special termination of STEP 2 of 
the algorithm was invoked; i.e., t* = 0 and the sample con- 
tained no zero values. For the case m^ = 5.0 and m 2 = 6.0, 
cycling between the following estimate values occurred. 




.0000 3.2191 9.6798 

.9310 2.9465 10.5754. 

The successive estimate pairs m^* and m 2 * produce values 
of t* which are on either side of the discontinuous ridge 
seen in Fig. 1 of section IV. As previously discussed, the 
efficiency function is bimodal in this region and causes the 
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above behavior. Since the asymptotic efficiency is high in 
this region, a larger sample size producing better initial 
estimates should resolve the inconsistency which results 
when cycling occurs. 
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VIII. TABLES 



Tables 1, 2 and 3 mentioned in the previous sections are 
contained on the following pages. Table 1 presents the op- 
timal t and efficiency values found for the range of m^ and 
m 2 from .01 - .09 (increments of .01), .1 - .9 (increments 
of .1), and 1.0 - 9.0 (increments of 1.0). Table 2 presents 
the optimal t and efficiency (relative to the efficiency of 
the Method of Moments) values found for the same range as 
used in Table 1. Table 3 presents the results of applying 
the algorithm of section V.B to the data sets of Appendix B. 
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TABLE is OPTIMAL T AND EFFICIENCY VALUES 
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EFF = BOTTCN VALUE OF EACH PAIR 



TABLE is OPTIMAL T AND EFFICIENCY VALUES 
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EFF = BOTTOM VALUE OF EACH PAIR 



TABLE l: OPTIMAL T AND EFFICIENCY VALUES 
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EFF = BOTTCP VALUE OF EACH PAIR 
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EFF » BOTTCM VALUE OF EACH PAIR 
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EFF = BOTTCM VALUE OF EACH PAIR 



TABLE is OPTIMAL T AND EFFICIENCY VALUES 
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EFF = BOTTOM VALUE OF EACH PAIR 
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TABLE i: OPTIMAL T AND EFFICIENCY VALUES 
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EFF = BOTTOM VALUE OF EACH PAIR 



TABLE l: OPTIMAL T AND EFFICIENCY VALUES 
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EFF = BOTTCM VALUE OF EACH PAIR 



TABLE 2: OPTIMAL T AND EFFICIENCY VALUES RELATIVE TO THE METHOD OF MOMENTS 
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EFF = BOTTOM VALUE OF EACH PAIR 



TABLE 3: NUMERICAL RESULTS 
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TABLE 3: NUMERICAL RESULTS 
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The sample contained no zero values and t* = 0, therefore the 
algorithm terminated with the Method of Moments estimators. 
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algorithm was terminated after 20 



TABLE 3: NUMERICAL RESULTS 
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TABLE 3: NUMERICAL RESULTS 
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IX. CONCLUSIONS 



The EPGF method provides estimators with asymptotic ef- 
ficiency close to that of the Maximum Likelihood estimators 
over most of the parameter values considered. The smallest 
value found was 54% (at m^ = 3.0 and m 2 = 9.0) and in most 
instances the asymptotic efficiency is greater than 97%. 

The EPGF method outperforms the Method of Moments in pro- 
ducing estimators of consistently higher asymptotic efficiency 
for all parameter values from .01 to 9.0. 

The efficiency of the EPGF method relative to the Method 
of Moments is less than 110% when m 2 < . 3 for most values of 
m^. As may be seen from Table 3, the relative efficiency 
drops below this value again for large m^ and m 2 beginning 
about m^ = 9.0 and m 2 = 3.0, and extending to m^ = 6.0 and 
m 2 = 9.0. 

In almost all cases, t* = 0. when m^ is less than 1.0. 

When this occurs the EPGF method reduces to estimation by the 
Method of Zero Frequency. The EPGF method becomes an exten- 
sion of the Zero Frequency method when the efficiency of the 
estimators of the latter decreases. This happens when t* > 

0. which is always the case for m^ > 3.0. 

Although the EPGF method is always better than the Method 
of Zero Frequency in terms of efficiency, it involves a 
greater computational effort when t* ^ 0. As shown in [4], 
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the Method of Moments estimators have higher asymptotic ef- 
ficiency than those of the Method of Zero Frequency whenever 
nu > 3.5 for all values of rr^ . The efficiency of the Method 
of Moments estimators is over 20% higher than the Zero Fre- 
quency estimators when m-^ > 4.5. Combined with the above 
discussion of the EPGF method, whenever t* > 0. the EPGF 
method will yield estimators of substantially higher asymp- 
totic efficiency than those found by the Zero Frequency method. 
The choice of which method to use always rests with the user. 
However, if it is assumed that m 2 > • 3 and that both m^ and 
m 2 are not large (greater than about 7.0) then the extra ef- 
fort required by the EPGF method is compensated for by a 

better than 10% increase in asymptotic efficiency of the 

« 

estimators. When it is supposed that m^ and m 2 are not to be 
found in the above region then estimation by the Method of 
Moments is recommended, as little increase in efficiency can 
be expected by using the EPGF method; i.e., going beyond 
STEP 1 of the algorithm presented in section V.B. 

As noted in section VII, the EPGF method as applied by 
the algorithm of section V.B can fail to converge. This is 
characterized by cycling between two different pairs of 
estimators when the sample size is small. When this occurs 
the remedy is to use a larger sample. The variance of a 
population having the Neyman Type A as an underlying distri- 
bution grows large as m^ and m 2 increase, especially when 
both are greater than 1.0. Assuming that the sample size is 
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large enough for the Central Limit theorem to provide a good 
approximation by a Normal distribution, the sample size re- 
quired to provide a reasonable confidence interval containing 
the true mean is quite large. Since x = m^rr^ is one of the 
estimating equations used in the EPGF method, it may be dif- 
ficult to provide good initial estimates for the iterative 
procedure when n is small. It has been found by experience, 
as the sample size is increased, the EPGF estimates tend to 
cluster more tightly than the Method of Moments estimators 
about the true parameter values. 

The EPGF method is not as easy to use as the Method of 
Moments, but with the aid of Table 2 and the use of the 
algorithm presented it can be applied readily with current 
programmable calculators. This is not always the case for 
estimation by Maximum Likelihood. 

Although the EPGF method appears to be time consuming 
and tedious, it is far easier to implement than the Maximum 
Likelihood method. By using Table 1, reasonably good esti- 
mators can be obtained with less effort than required by the 
latter method. 
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APPENDIX A 



COMPUTER PROGRAMS 



Four computer programs were used in the compilation of 

this paper and are given on the following pages. All were 

programmed in the standard Fortran IV language and in double 

precision. Program 1 was used to compute the values of the 

information determinant | A | of section III for m^ and m 2 in 

the range of Table 1. Program 2 was run to compute the 

values of t* and EFF(t*) and produce Table 1. Table 2 was 

produced by slightly modifying program 2 to compute and table 

the efficiency ratio EFF/EFF of section VI. Program 3 was 

mm 

used to perform the algorithm of section V.B and produce the 
estimates for the numerical examples of section VII, which 
are found in Table 3. Program 4 was used in simulating ran- 
dom numbers from the Neyman Type A distribution. The sub- 
routine LP0IS1 found in program 4 is a local random number 
generator for Poisson distributed random deviates. 

Programs 1, 2, and 3 are complete and may be used di- 
rectly. Program 4 must be modified to be compatible with 
the user's local random number generation subroutines. 
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non non 



COMPUTER PROGRAM 1 



C COMPUTATION OF INFORMATION DETERMINANT GIVEN Ml AND M2 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL *8 Ml, M2 
COMMON Ml,M2,BETA,DET 
DIMENSION PARM(27) , DETQ (27,27) 

DATA PARM/1 . D- 2 , 2 . D- 2 , 3 . D- 2 , 4 . D- 2 , 5 . D- 2 , 6 . D-2 , 7 . D- 2 , 

* 8 . D-2 , 9 . D-2 , 1 . D-l , 2 . D-l , 3 . D-l , 4 . D-l , 5 . D-l , 6 . D-l , 7 . D-l , 

* 8 . D- 1 , 9 . D-l , 1 . DO , 2 . DO , 3 . DO , 4 . DO , 5 . DO , 6 . DO , 7 . DO , 8 . DO , 

* 9. DO/ 

CREATE A TABLE OF DETERMINANT VALUES FOR VARYING 
PARAMETERS OF 
Ml AND M2 
J1 = 27 
J2 = 27 

DO 999 I = 27,27 
Ml = PARM(I) 

DO 99 J = J1,J2 
M2 = PARM(J) 

BETA = Ml *DEXP (-M2 ) 

VALUES HAVE BEEN SET FOR Ml AND M2 

COMPUTE INFORMATION DETERMINANT 
CALL DETQS 
DETQ ( I , J) = DET 
99 CONTINUE 

IF (I.NE.l) GO TO 150 
WRITE (6,101) (PARM(L) ,L=J1, J2) 

101 FORMAT (5X,3 (3X,F12.10) ) 

150 WRITE (6,151) Ml , (DETQ ( I , L) , L=J1 , J2 ) 

151 FORMAT (F5.2,3(3X,F12.10) ) 

999 CONTINUE 

STOP 

END 

C 

SUBROUTINE DETQS 
IMPLICIT REAL *8 (A-H,0-Z) 

REAL* 8 Ml, M2 
COMMON Ml, M2, BETA, DET 
C COMPUTE Q GIVEN Ml AND M2 

H = l.D-7 

PO = DEXP (BETA - Ml) 

RX = BETA 
QX = BETA*RX 
Q = QX 

ELN = DLOG (BETA) + BETA 
DO 10 IX = 1,500 
X = DFLOAT ( IX) 

ELNEXT = SERIES (X+l. DO) 
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RNEXT = ESXP (ELNEXT - ELN) 

ELN = ELNEXT 

QX = QX*M2* (RNEXT/X) * (RNEXT/RX) 

RX = RNEXT 

IF (QX/Q .LT. H) GO TO 20 
K = IX 

10 Q = Q + QX 

20 Q = M2*M2*PO*Q 

A = M1*M2*M2*M2 
B = M1*M2*M2* (Ml+M2+M1*M2) 

DET = (1.D0 + M2 ) *Q/A - B/A 
RETURN 
END 
C 

FUNCTION SERIES (X) 

IMPLICIT REAL*8 (A-H, O-Z) 

REAL* 8 Ml, M2 

COMMON Ml, M2, BETA, DET 

SERIES = 0 . DO 

C FIND INDEX OF THE MAX TERM: AIM 

AIM = 1.D0 
Cl = 1 . DO 
C2 = 2. DO 
DO 10 K = 1,100 

FNUM = C2*AIM*DLOG (BETA/AIM) + C2*X - Cl 

FDEN = C2 * (X + AIM) - Cl 

AIN EXT = AIM* (Cl + FNUM/FDEN) 

IF (AINEXT .LE. l.D-5) GO TO 11 
AIM = AINEXT 

10 CONTINUE 

11 IF (AIM .LE. 1 . DO ) AIM = 1 . DO 

IF (AIM .GT. 1 . DO ) AIM = DFLOAT ( IDINT (AIM + .5D0)) 

IM = IDINT (AIM) 

C COMPUTE MAX TERM OF THE SERIES 

ALNIM = 0 . DO 
DO 20 I = 1 , IM 
AI = DFLOAT (I) 

20 ALNIM = ALNIM + DLOG(AI) 

ALNIM = X*DLOG (AIM) + AIM*DLOG (BETA) - ALNIM 
C COMPUTE EACH TERM (K) OF SERIES AND SUM TO THE SERIES 

DO 40 K = 1,300 
AK = DFLOAT (K) 

ALNK = 0.D0 
DO 30 J = 1 , K 
AJ = DFLOAT (J) 

30 ALNK = ALNK + DLOG(AJ) 

ALNK = X*DLOG (AK) + AK*DLOG (BETA) - ALNK 
BK = ALNK - ALNIM 

IF (BK .GE. -161. 1D0) SERIES = SERIES + DEXP (BK) 

IF (BK .LT. -161. 1D0 .AND. K . GT . IM) GO TO 50 
40 CONTINUE 

C COMPUTE THE NATURAL LOG OF THE XTH MOMENT. 

50 SERIES = ALNIM + DLOG (SERIES) 

RETURN 

END 
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COMPUTER PROGRAM 2 



C 



C 

C 

C 

C 

C 

C 

C 



50 

101 

99 

C 

C 

C 



199 

C 



297 



COMPUTATION OF OPTIMAL T AND EFFICIENCY GIVEN Ml AND M2 
IMPLICIT REAL *8 (A-H,0-Z) 

REAL *8 Ml, M2 

COMMON Ml , M2 , DET , TSTAR , ESTAR 

DIMENSION PARM (27) ,TOPT(27,27) ,EOPT(27,27) ,DETQ(27,27) 
CREATE A TABLE OF T-OPTIMAL VALUES FOR VARYING 
PARAMETERS OF 
Ml AND M2 



READ IN Ml AND M2 AND CORRESPONDING INFORMATION 
DETERMINANT VALUES FROM A DISK 



J2 = 0 

DO 99 K = 1,9 

J1 = J2 + 1 

J2 = J1 + 2 

DO 99 I = 1,27 

IF ( I .NE. 1) GO TO 50 

READ (5,101) (PARM(J) ,J=J1,J2) 

READ (5,101) (DETQ(I,J) ,J=J1,J2) 
FORMAT (5X, 3 (3X,F12 . 10) ) 

CONTINUE 

SET Ml AND M2 VALUES FROM PARM 
COMPUTE OPTIMAL T-STAR = TOPT AND 
E-STAR = EOPT 
DO 199 I = 1,27 
Ml = PARM (I) 

DO 199 J = 1,27 
M2 = PARM(J) 

DET = DETQ ( I , J) 

CALL OPTIML 
TOPT (I , J) = TSTAR 
EOPT (I , J) = ESTAR 
CONTINUE 
PRINT OUT 
IU 
DO 
IL 
IU 
JU 
DO 
JL 
JU 



OPTIMAL EFFICIENCY 



= 0 

299 KI 
= IU + 
= IL + 
= 0 

298 KJ 
= JU + 
= JL + 



TABLE OF OPTIMAL T AND EFFICIENCY 

= 1,3 
1 
8 

= 1,3 
1 
8 



WRITE (6,201) 

WRITE (6,202) 

WRITE (6,203) 

WRITE (6,204) (PARM(J) , J=JL,JU) 

DO 297 I = IL , IU 

WRITE (6,205) PARM (I ) , (TOPT ( I , J) , J=JL, JU) 
WRITE (6,206) (EOPT (I , J) , J=JL, JU) 

CONTINUE 



66 



WRITE (6,207) 

WRITE (6,208) 

WRITE (6,209) 

298 CONTINUE 

299 CONTINUE 

201 FORMAT (' 1' ,//////) 

202 FORMAT ( ' O' , 33X, 'TABLE 1: OPTIMAL T AND EFFICIENCY 

* VALUES',//) 

203 FORMAT ( ' 0 ' , 9X , ' Ml VALUES ', 35X, ' M2 VALUES') 

204 FORMAT ('O' ,19X,9(2X,F6.4) ) 

205 FORMAT ( 'O' , 11X,F4 . 2 , 4X, 9 (2X,F6 . 4) ) 

206 FORMAT ( ' ' , 19X, 9 (2X,F6 . 4) ) 

207 FORMAT (' O' ,/,22X, 'FOR EACH ENTRY IN THE TABLE:') 

208 FORMAT ( 'O' ,21X, ' T* = TOP VALUE OF EACH PAIR') 

209 FORMAT ('O' , 2 IX , ' EFF = BOTTOM VALUE OF EACH PAIR') 
STOP 

END 

C 

SUBROUTINE OPTIML 
IMPLICIT REAL * 8 ( A- H , O- Z ) 

REAL* 8 Ml, M2 

COMMON Ml , M2 , DET , TSTAR , ESTAR 
FLAG = 0 

THETA = 9999. D-4 
E2 = 0 . DO 

C FIND A COARSE BRACKET AROUND THETA: (A,B) . 

5 THETA = THETA - l.D-2 

IF (THETA .LT. 0 . DO ) GO TO 15 
ET = EFF (THETA) 

IF (ET .LT. E2 ) GO TO 15 
E2 = ET 
GO TO 5 

15 A = THETA 

IF (A .LT. 0 . DO ) A = 0 . DO 
B = THETA + 2.D-2 
IF (B .GT. 9999. D-4) B = 9999. D-4 
C 

PERFORM A GOLDEN SECTION SEARCH FOR TSTAR TO MAX 
C EFF (TSTAR) 

C WITH STARTING INTERVAL (A, B) . 

R = (-1.D0 + DSQRT ( 5 . DO ) ) /2 . DO 
20 Al = ( 1 . DO - R) * (B - A) + A 
B1 = R* (B - A) + A 
EA1 = EFF (Al) 

EBl = EFF(Bl) 

IF (DABS (EAl - EBl) .LT. l.D-7) GO TO 30 
IF (EAl .GT. EBl) B = Bl 
IF (EAl .LT. EBl) A = Al 
GO TO 20 

30 TSTAR = (Al + B1)/2.D0 
ESTAR = EBl 

C CHECK FOR A MAX FROM THE LEFT END 

IF (FLAG .EQ. 1) GO TO 40 
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IF (TSTAR .LT. l.D-1) RETURN 
T = TSTAR 
E = ESTAR 
THETA = O.DO 
E2 = O.DO 

35 THETA = THETA + l.D-2 

IF (THETA .GT. 99.D-2) RETURN 
ET = EFF (THETA) 

IF (ET .LT. E2 ) GO TO 36 
E2 = ET 
GO TO 35 

36 A = THETA - 2.D-2 

IF (A .LT. O.DO) A = O.DO 
B = THETA 
FLAG = 1 
GO TO 20 

40 IF (E .LT. ESTAR) RETURN 
TSTAR = T 
ESTAR = E 
RETURN 
END 
C 

FUNCTION EFF (T) 

IMPLICIT REAL* 8 ( A-H , 0- Z ) 

REAL* 8 Ml, M2 

COMMON Ml, M2, DET, TSTAR, ESTAR 
C COMPUTE PROBABILITY GENERATING FUNCTION G (T) 

F = DEXP (M2*T - M2) 

G = DEXP (-Ml* (1. DO - F) ) 

C COMPUTE G (T* *2 ) 

F2 = DEXP (M2 *T*T - M2) 

G2 = DEXP (-Ml* (1. DO - F2)) 

C COMPUTE NUMERATOR OF EFF (T) 

FN = 1 . DO + (M2*T - M2 - l.D0)*F 
FN = Ml*G*G*FN*FN 

C COMPUTE DENOMINATOR OF EFF (T) 

FD = 1 . DO + M2 + M1*M2*(T*F - l.D0)*(T*F - 1.D0) 
FD = M2* ( (1. DO + M2) *G2 - G*G*FD) 

FD = DET*FD 

C COMPUTE EFFICIENCY EFF (T) 

EFF = FN/FD 

RETURN 

END 
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COMPUTER PROGRAM 3 



C ESTIMATION OF Ml AND M2 

C 

IMPLICIT REAL* 8 (A-H,0-Z) 

DIMENSION IX (450, 10) , M2 (9) 

COMMON HATM1 , HATM2 
READ (5,9) (M2 (J) ,J=1,9) 

9 FORMAT (3X, 915) 

C READ IN MATRIX OF RANDOM NOS . : IX 

READ (5,10) ( (IX(I,J) ,J=1, 10) ,1=1,450) 

10 FORMAT (13, 915) 

NNN = 50 

ANN = DFLOAT (NNN) 

IFLAG = -1 
DO 90 I = 1,9 
IFLAG = IFLAG* (-1) 

IK = NNN* (1-1) + 1 
IKE = IK + NNN - 1 
PM1 = DFLOAT (IX (IK, 1) ) 

IF (IFLAG .LT. 0) GO TO 15 
WRITE (6,100) 

100 FORMAT (' 1' ,////, T37, 'TABLE 3: NUMERICAL RESULTS') 

WRITE (6,101) 

101 FORMAT ('O' , / , T19 , ' Ml ' , T2 7 , ' M2 ' , T35 , ' Ml ' , T4 3 , ' M2 ' , T51 , 

* ' T* ' , T58 , ' T* 1 , T65 , ’Ml*’ ,T73, 'M2*' ,T80, 'NO. * ) 

WRITE (6,102) 

102 FORMAT ( ' ' , T34 , ' (MM) ' , T42 , ' (MM) ' , T50 , ' INIT ' , T56 , ' FINAL ' , 

* T80 , ' ITR. ' ) 

15 CONTINUE 

DO 90 J = 2,10 

JK = J - 1 

PM2 = DFLOAT (M2 (JK) ) 

XBAR = 0.D0 

XSQ = 0.D0 

DO 20 K = IK, IKE 

XBAR = XBAR + DFLOAT (IX (K, J) ) 

20 XSQ = XSQ + DFLOAT (IX (K,J) *IX(K,J) ) 

XBAR = XBAR/ ANN 

XVAR = (XSQ - ANN* XBAR* XBAR) / (ANN - 1 . DO ) 

HATM2 = (XVAR - XBAR) /XBAR 

HATM1 = XBAR/HATM2 

AMI = PM1 

AM2 = PM2 

BMl = HATMl 

BM2 = HATM2 

C FIND ESTIMATORS USING EMPIRICAL PGF 

PAST = HATM2 
JJJ = 0 

30 CALL TOPT (T) 

IF (JJJ .EQ. 0) TI = T 
JJJ = JJJ + 1 
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IF (JJJ .GT. 20) GO TO 65 
TBAR = 0.D0 
DO 60 KK = IK, IKE 

IF (IX (KK , J) .EQ. 0) TBAR = TBAR + 1 . DO 

IF (T .GT. l.D-4 .AND. IX(KK,J) .NE. 0) TBAR = TBAR + 

* T* *IX (KK , J) 

60 CONTINUE 

IF (TBAR .LE. l.D-4) GO TO 70 
TBAR = TBAR/ ANN 

CALL ANEWT (TBAR, PAST , XBAR, T , PM2 ) 

HATMl = XBAR/PM2 

IF (DABS (HATM2 - PM2) . LE . l.D-4) GO TO 80 
HATM2 = PM2 
GO TO 30 

65 WRITE (6,103) AMI , AM2 , BMl , BM2 , TI 

103 FORMAT ('O' , T14 , 4F8 . 3 , T49 ,F5 . 4 ,T56 , ' ' ,T64, ' 

* T72, ' ' , T81 , ' + ' ) 

GO TO 90 

70 WRITE (6,104) AMI , AM2 , BMl , BM2 , TI , T , BMl , BM2 , JJJ 

104 FORMAT ('O' , T14 , 4F8 . 3 , T49 ,F5 . 4 , T56 , F5 . 4 , 2F8 . 3 , 15 , '*') 
GO TO 90 

80 WRITE (6,105) AMI ,AM2 , BMl , BM2 , TI , T , HATMl , PM2 , JJJ 

105 FORMAT ( 'O' , T14 , 4F8 . 3 , T49 , F5 . 4 , T56 , F5 . 4 , 2F8 . 3 , 15 ) 

90 CONTINUE 

WRITE (6,106) 

106 FORMAT ('1') 

STOP 

END 

C 

SUBROUTINE ANEWT (TBAR, YIN , XBAR , T , YOUT) 

C 

IMPLICIT REAL* 8 (A-H,0-Z) 

COMMON HATMl, HATM2 
TL = DLOG(TBAR) 

DO 20 KKK = 1,20 

GY = DEXP (YIN* (T-1.D0) ) 

FY = YIN*TL + XBAR* (1. DO - GY) 

FPY = TL - XBAR* GY* (T - 1.D0) 

YOUT = YIN - FY/FPY 

IF (DABS (YOUT - YIN) .LE. l.D-4) GO TO 99 
YIN = YOUT 
20 CONTINUE 

99 RETURN 

END 
C 

SUBROUTINE TOPT (T) 

C 

IMPLICIT REAL * 8 ( A- H , O- Z ) 

COMMON HATMl, HATM2 
FLAG = 0 

THETA = 9999. D-4 
E2 = 0 . DO 

C FIND A COARSE BRACKET AROUND THETA: (A,B) 
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5 THETA = THETA - l.D-2 

IF (THETA .LT. O.DO) GO TO 15 
ET = EFF (THETA) 

IF (ET .LT. E2 ) GO TO 15 
E2 = ET 
GO TO 5 

15 A = THETA 

IF (A .LT. O.DO) A = O.DO 
B = THETA + 2.D-2 
IF (B .GT. 9999. D-4) B = 9999. D-4 
R = (-1.D0 + DSQRT ( 5 . DO ) ) /2 . DO 
20 A1 = (1 . DO - R) * (B - A) + A 
B1 = R* (B - A) + A 
EA1 = EFF(Al) 

EB1 = EFF(Bl) 

IF (DABS (EA1-EB1) . LT . l.D-15) GO TO 30 

IF (EA1 .GT. EB1) B = B1 
IF ( EAl .LT. EB1) A = A1 
GO TO 20 

30 T = ( A1 + B1)/2.D0 
C CHECK FOR MAX FROM THE LEFT END 

IF (FLAG .EQ. 1) GO TO 40 
IF ( T .LT. l.D-1) RETURN 
TSTAR = T 
ESTAR = EB1 
THETA = O.DO 
E2 = O.DO 

35 THETA = THETA + l.D-2 

IF (THETA .GT. 99.D-2) RETURN 
ET = EFF (THETA) 

IF (ET .LT. E2 ) GO TO 36 
E2 = ET 
GO TO 35 

36 A = THETA - 2.D-2 

IF (A .LT. O.DO) A = O.DO 
B = THETA 
FLAG = 1 
GO TO 20 

40 IF (ESTAR .LT.EB1) RETURN 

T = TSTAR 
RETURN 
END 
C 

FUNCTION EFF (Q) 

C 

IMPLICIT REAL * 8 ( A- H , O- Z ) 

COMMON HATM1 , HATM2 
AMI = HATMl 
AM2 = HATM2 
F = DEXP ( AM2 *Q - AM2 ) 

G = DEXP (-AM1* (1.D0 - F) ) 

F2 = DEXP (AM2*Q*Q - AM2 ) 

G2 = DEXP (-AM1* (1.D0 - F2) ) 
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FN = 1 . DO + (AM2*Q - AM2 - l.DO)*F 
FN = AMl*G*G*FN*FN 

FD = 1 . DO + AM2 + AM1*AM2*(Q*F - l.DO)*(Q*F - 1 . DO ) 
FD = AM2* ( (l.DO + AM2)*G2 - G*G*FD) 

EFF = FN/FD 

RETURN 

END 
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COMPUTER PROGRAM 4 



C 

C 



1 

2 

3 

10 

C 

99 

50 

53 

100 
55 
101 
200 



RANDOM NUMBER GENERATION FROM NEYMAN ' S TYPE A 
DISTRIBUTION 
DIMENSION AD (50) 

DATA ISEED1,ISEED2,IS0RT,KSIZE/431219, 187648 ,0,50/ 
CALL OVFLOW 
DO 10 Ml = 1,9 
PMl = FLOAT (Ml) 

CALL LPOIS1 ( ISEED1 , AD , KSI ZE , PMl , I SORT) 

WRITE (6,1) 

FORMAT Cl', ///////) 

WRITE (6,2) 

FORMAT ( ' ' ,14X, ' RANDOM NOS. FROM THE NEYMAN TYPE A 

* DISTRIBUTION') 

WRITE (6,3) 

FORMAT ('O' , // , 12X ' Ml VALUE ', T40 ,' M2 VALUES') 

CALL DGHTER(AD,ISEED2,M1,KSIZE) 

CONTINUE 

STOP 

END 



SUBROUTINE DGHTER (AD , ISEED2 ,Ml) 
DIMENSION AD (50) , Y (100) ,IX(50,9) 
WRITE (6,99) ( J , J=1 , 9 ) 

FORMAT ('O' , 18X, 915 ,//) 

DO 200 K = 1,50 
ISIZE = INT (AD (K) ) 

IF (ISIZE .NE. 0) GO TO 53 

DO 50 M2 = 1,9 

IX (K,M2 ) = 0 

GO TO 55 

DO 100 M2 = 1,9 

PM2 = FLOAT (M2) 

CALL LPOISl ( I SEED 2 , Y , ISI ZE , PM2 , 0 ) 

IX (K ,M2 ) = 0 

DO 100 I = 1, ISIZE 

IX (K , M2 ) = IX (K ,M2 ) + INT (Y (I) ) 

WRITE (6,101) Ml, (IX(K,J) ,J=1,9) 

FORMAT ( ' ' ,10X,I5,3X,9I5) 

CONTINUE 

RETURN 

END 
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APPENDIX B 



RANDOM NUMBERS FROM THE 
NEYMAN TYPE A DISTRIBUTION 



Computer program 4 was used to generate the eighty-one 
sets of random numbers for the numerical examples of Table 3. 
The parameter range used was m^ and m2 from 1 to 9 in incre- 
ments of 1 and the sample size was 50 in all cases. The 
following pages contain the data sets generated. 
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RANDOM NOS 



FROM THE NEY1 AN TY D E A DISTRIBUTION 



Ml VALUE 



2 



3 



M2 VALUES 

A 5 6 7 8 9 



1 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

i 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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RANDOM NCS 



FROM THE NEYMAN TYPE A DISTRIBUTION 



Ml VALUE 

12 3 



COO 

c 0 c 

4 5 6 

0 0 0 

4 5 7 

0 0 0 

3 3 3 




2 5 8 

5 5 6 

6 5 15 

0 8 3 

0 C 0 

2 5 8 

0 0 0 

5 9 25 

2 4 6 

2 4 4 

C 2 2 



1 3 4 
1 2 5 
0 0 0 
0 5 3 
7 5 9 
0 0 0 



VALUES 



5 


6 


7 


8 


9 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


21 


14 


26 


21 


25 


0 


0 


0 


0 


0 


25 


21 


34 


35 


33 


0 


0 


0 


0 


0 


13 


6 


16 


15 


15 


10 


18 


10 


19 


17 


20 


27 


28 


19 


29 


8 


14 


11 


11 


23 


11 


16 


13 


20 


18 


19 


22 


21 


27 


33 


30 


24 


20 


31 


45 


0 


0 


0 


0 


0 


14 


17 


13 


21 


34 


0 


0 


0 


0 


0 


7 


14 


1 1 


19 


10 


17 


23 


24 


20 


29 


32 


28 


36 


43 


49 


0 


0 


0 


0 


0 


12 


9 


22 


12 


22 


7 


6 


4 


3 


10 


15 


26 


34 


39 


34 


8 


15 


1 1 


14 


17 


9 


7 


8 


9 


7 


20 


20 


24 


35 


35 


14 


19 


12 


18 


14 


2 


5 


4 


8 


8 


13 


11 


6 


14 


22 


18 


11 


28 


16 


23 


10 


7 


14 


13 


18 


13 


26 


24 


28 


36 


12 


18 


13 


17 


15 


0 


0 


0 


0 


0 


10 


9 


1 5 


13 


18 


0 


0 


0 


0 


0 


40 


39 


51 


67 


67 


9 


15 


11 


16 


22 


8 


4 


4 


10 


10 


4 


9 


5 


12 


13 


0 


0 


0 


0 


0 


13 


19 


21 


25 


26 


9 


14 


13 


14 


22 


1 


6 


9 


1C 


9 


9 


19 


11 


10 


1 5 


6 


5 


7 


8 


7 


0 


0 


0 


0 


0 


8 


5 


9 


4 


12 


12 


15 


19 


28 


44 


0 


0 


0 


C 


0 



M2 

4 

0 

0 

17 

0 

13 

0 

2 

10 

18 

10 

6 

20 

20 

0 

7 

0 

13 

11 

13 

0 

12 

3 

17 

6 

4 

13 

7 

1 

10 

10 

5 

14 

10 

0 

7 

0 

31 

14 

6 

3 

0 

18 

9 

5 

8 

l 

0 

3 

14 

0 
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RANDOM NOS. FROM THE NEYMAN TYPE A DISTRIBUTION 



Ml VALUE 

12 3 



4 6 14 

4 14 13 

4 1C 10 

4 3 5 

6 7 17 

15 5 

12 7 

3 7 7 

3 14 11 

0 0 4 

3 7 10 

COO 

1 6 2 

3 14 25 

7 9 17 

3 5 7 

3 7 13 

3 7 8 

16 8 

8 12 20 

4 4 12 

2 9 15 

2 6 11 




C 5 6 
0 0 0 
3 7 6 
1 6 11 



VALUES 



5 


6 


7 


8 


9 


22 


37 


29 


30 


31 


24 


40 


30 


44 


49 


19 


21 


20 


lo 


25 


1 1 


10 


9 


20 


16 


25 


43 


44 


58 


63 


8 


12 


14 


16 


22 


1 7 


12 


24 


24 


28 


14 


9 


12 


15 


15 


22 


36 


48 


54 


64 


0 


4 


7 


7 


9 


22 


20 


34 


36 


40 


0 


0 


0 


0 


0 


5 


5 


6 


10 


9 


32 


42 


49 


54 


41 


31 


33 


38 


40 


46 


15 


11 


21 


24 


26 


24 


25 


28 


34 


35 


1 1 


16 


1 5 


30 


27 


16 


19 


17 


29 


23 


23 


26 


33 


26 


26 


23 


3 L 


32 


22 


40 


16 


19 


23 


22 


21 


15 


15 


3 1 


16 


30 


4 


5 


1 


10 


8 


15 


19 


1 6 


25 


24 


23 


33 


30 


32 


33 


0 


6 


3 


8 


10 


6 


10 


8 


8 


10 


18 


25 


31 


23 


38 


6 


9 


9 


8 


7 


7 


13 


16 


18 


22 


5 


8 


4 


6 


14 


12 


7 


15 


16 


22 


16 


22 


18 


20 


28 


13 


13 


19 


25 


26 


20 


21 


25 


27 


25 


21 


30 


32 


42 


46 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


9 


8 


7 


1C 


16 


0 


0 


0 


0 


0 


25 


19 


14 


26 


20 


18 


19 


26 


36 


41 


16 


27 


26 


29 


44 


25 


27 


43 


37 


40 


30 


38 


39 


50 


49 


12 


18 


14 


16 


18 


0 


0 


0 


0 


0 


23 


27 


27 


27 


35 


25 


20 


26 


31 


41 



M2 

4 

19 

21 

10 

9 

29 

9 

19 

5 

20 

7 

10 

0 

5 

19 

18 

14 

18 

7 

16 

16 

15 

11 

14 

4 

14 

20 

3 

5 

18 

1 

io 

2 

14 

12 

16 

10 

12 

0 

0 

8 

0 

11 

14 

20 

27 

32 

8 

0 

16 

25 



77 



4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

A. 



.NOS. FROM THE NEV MAN TYPE A DISTRIBUTION 



1 



2 

3 

7 
1 

4 
3 
3 

8 
1 
0 
2 
9 
7 

3 
13 

7 

5 

4 
3 
2 

5 
1 
1 

5 

8 
2 

2 

£ 

2 

2 

2 

3 

2 

•3 

2 

2 

0 

2 

I 

1 

2 
3 
0 

0 

a 

w 

6 
9 

1 2 
5 
0 



M2 VALUES 



2 


3 


4 


5 


6 


7 


8 


9 


8 


1 6 


29 


17 


24 


31 


29 


44 


8 


15 


22 


23 


20 


1 7 


28 


38 


12 


23 


30 


30 


40 


37 


49 


39 


8 


7 


7 


10 


13 


15 


15 


1 5 


8 


10 


26 


22 


28 


26 


32 


39 


10 


17 


17 


34 


33 


25 


51 


76 


7 


1 9 


15 


32 


40 


32 


45 


60 


17 


26 


25 


30 


41 


73 


63 


73 


5 


8 


10 


14 


19 


19 


20 


26 


13 


5 


11 


16 


23 


22 


20 


33 


2 


6 


4 


3 


6 


7 


10 


12 


10 


14 


17 


24 


31 


22 


26 


44 


10 


1 4 


26 


18 


26 


31 


44 


49 


5 


6 


14 


11 


18 


1 6 


21 


22 


9 


26 


22 


20 


42 


47 


69 


58 


6 


12 


21 


19 


33 


44 


41 


49 


15 


21 


29 


27 


30 


50 


58 


52 


4 


10 


3 


10 


8 


12 


14 


19 


9 


6 


10 


1C 


15 


1 1 


25 


26 


9 


12 


14 


21 


29 


39 


30 


35 


20 


20 


36 


48 


54 


66 


69 


83 


9 


7 


13 


1 1 


28 


27 


16 


22 


2 


1 


5 


7 


2 


3 


9 


9 


1 


12 


12 


13 


11 


16 


14 


15 


10 


3 1 


26 


52 


44 


56 


76 


65 


4 


5 


5 


11 


7 


10 


19 


1 6 


1 


7 


7 


12 


14 


12 


10 


17 


9 


13 


20 


13 


25 


33 


43 


37 


6 


13 


17 


12 


23 


27 


37 


30 


6 


9 


13 


10 


11 


13 


19 


17 


7 


13 


14 


1 6 


31 


2 1 


31 


41 


1 


7 


11 


14 


10 


24 


23 


25 


4 


12 


12 


12 


17 


19 


23 


25 


8 


2 C 


13 


20 


41 


29 


41 


30 


6 


5 


6 


10 


12 


10 


14 


13 


3 


6 


6 


9 


11 


17 


13 


15 


5 


14 


8 


10 


10 


1 8 


21 


10 


7 


6 


10 


15 


22 


1 1 


27 


27 


14 


20 


21 


28 


55 


52 


61 


55 


2 


6 


12 


14 


10 


22 


7 


19 


7 


7 

4 


11 


12 


12 


20 


27 


26 


15 


16 


23 


32 


38 


52 


41 


54 


7 


6 


7 


10 


18 


20 


17 


14 


1 


1 


3 


6 


7 


6 


3 


10 


4 


8 


10 


1 0 


6 


8 


27 


10 


12 


22 


25 


38 


37 


37 


48 


59 


12 


17 


13 


29 


39 


23 


47 


59 


14 


34 


20 


36 


53 


44 


69 


3 6 


6 


10 


24 


17 


21 


36 


34 


39 


2 


4 


13 


9 


11 


14 


20 


21 



78 



in u > in m u y m in u > >n in m m u 'in in in in in u ' in m m in tn in in u in in in in in *n in in u 1 in in in m u > in tn 1 u ' m in « ' to tn in 



RANDOM NOS. FROM 



HE NEYMAN TYPE A DISTRIBUTION 



Ml VALUE 



1 


2 


3 


4 


17 


20 


6 


14 


1 5 


9 


18 


15 


1 1 


14 


27 


6 


11 


29 


11 


10 


23 


4 


6 


7 


C 


0 


0 


7 


8 


22 


9 


24 


23 


2 


9 


12 


8 


10 


16 


3 


3 


2 


3 


6 


5 


0 


3 


2 


3 


13 


15 


•a 


12 


13 


3 


8 


13 


1 C 


23 


33 


2 


0 


2 


6 


5 


16 


c 


10 


16 


5 


10 


12 


7 


1C 


16 


2 


5 


13 


2 


7 


20 


7 


16 


30 


2 


5 


11 


5 


14 


18 


8 


12 


2 9 


1 


8 


4 


4 


29 


24 


4 


13 


1 7 


4 


13 


16 


6 


4 


12 


14 


19 


3 1 


5 


14 


16 


7 


6 


19 


5 


11 


13 


12 


20 


34 


4 


5 


15 


0 


3 


3 


4 


16 


13 


8 


20 


24 


4 


2 


4 


8 


20 


29 


1 


4 


11 


1 


8 


8 


5 


1 


12 


9 


17 


30 



VALUES 



5 


6 


7 


8 


9 


36 


37 


39 


46 


70 


39 


36 


51 


55 


50 


51 


34 


41 


62 


72 


43 


46 


49 


50 


60 


3 0 


37 


36 


41 


66 


29 


28 


37 


46 


56 


11 


10 


21 


25 


28 


0 


0 


0 


0 


0 


33 


34 


46 


C -3 
-s — ■ 


47 


40 


63 


53 


73 


80 


12 


19 


24 


25 


27 


29 


35 


49 


58 


60 


9 


7 


5 


8 


9 


4 


8 


1 5 


14 


16 


10 


0 


5 


6 


7 


27 


28 


36 


32 


47 


29 


39 


44 


62 


64 


23 


40 


29 


40 


44 


61 


49 


69 


73 


97 


6 


2 


12 


17 


10 


21 


34 


42 


39 


34 


25 


31 


32 


41 


35 


29 


27 


42 


47 


51 


32 


48 


43 


47 


67 


18 


17 


23 


26 


25 


29 


25 


41 


32 


47 


42 


60 


48 


78 


73 


22 


29 


2 6 


36 


41 


24 


38 


52 


67 


53 


5 1 


52 


58 


79 


33 


14 


5 


16 


19 


17 


45 


44 


62 


75 


95 


34 


43 


56 


48 


69 


29 


24 


37 


43 


47 


17 


17 


25 


29 


30 


26 


41 


49 


63 


84 


22 


32 


35 


45 


46 


39 


54 


59 


66 


66 


26 


36 


50 


50 


57 


39 


58 


59 


74 


77 


20 


27 


24 


35 


30 


6 


3 


5 


7 


11 


30 


26 


37 


37 


41 


29 


51 


51 


65 


74 


3 


15 


12 


17 


15 


50 


58 


87 


88 


82 


6 


6 


9 


18 


10 


13 


18 


17 


24 


22 


27 


32 


49 


35 


43 


45 


55 


6 3 


100 


96 



M2 

4 

22 

29 

26 

34 

39 

23 

13 

0 

20 

33 

15 

35 

4 

9 

7 

22 

27 

17 

50 

8 

17 

18 

16 

17 

14 

22 

25 

17 

34 

37 

10 

26 

23 

20 

22 

23 

22 

33 

28 

34 

21 

4 

19 

25 

4 

48 

10 

11 

24 

42 



79 



6 

£ 

6 

6 

£ 

6 

6 

£ 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

£ 

6 

6 

£ 

£ 

6 

6 

6 

6 

£ 

6 

6 

£ 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

£ 

6 

6 

6 

6 

6 

6 

6 



NCS . FROM ThE NEYMAN TYPE A DISTRIBUTION 

^2 VALUES 



1 


2 




4 


5 


6 


7 


8 


9 


5 


17 


16 


24 


46 


29 


41 


55 


55 


2 


5 


7 


5 


7 


15 


17 


13 


14 


5 


5 


13 


15 


22 


22 


19 


20 


17 


6 


12 


19 


23 


25 


46 


47 


67 


66 


3 


6 


16 


10 


18 


28 


26 


39 


41 


6 


10 


6 


19 


32 


21 


37 


39 


51 


7 


16 


1 1 


24 


33 


30 


36 


50 


54 


2 


5 


15 


15 


30 


21 


30 


34 


27 


4 


18 


20 


29 


40 


42 


66 


68 


73 


5 


9 


a 


17 


17 


20 


21 


21 


3 1 


0 


3 


2 


12 


9 


14 


12 


22 


17 


4 


9 


24 


29 


40 


52 


47 


58 


64 


4 


9 


13 


12 


20 


21 


25 


25 


33 


7 


11 


24 


36 


33 


52 


43 


66 


6 1 


3 


11 


16 


22 


33 


28 


37 


30 


49 


2 


13 


16 


20 


26 


36 


31 


45 


57 


2 


9 


3 


8 


16 


26 


25 


28 


27 


5 


13 


11 


21 


20 


29 


36 


40 


49 


1C 


24 


32 


39 


49 


49 


72 


67 


76 


13 


22 


29 


27 


54 


56 


79 


78 


102 


4 


11 


22 


24 


18 


30 


32 


39 


41 


2 


17 


15 


21 


33 


33 


33 


46 


36 


4 


6 


14 


18 


19 


23 


30 


30 


32 


1 


14 


19 


27 


25 


39 


49 


41 


57 


3 


5 


15 


25 


20 


32 


28 


36 


46 


7 


7 


9 


20 


39 


29 


32 


46 


7 6 


7 


10 


1 1 


21 


33 


28 


34 


46 


44 


7 


10 


1 8 


26 


37 


31 


45 


43 


65 


5 


11 


23 


30 


34 


46 


53 


52 


61 


1 6 


15 


34 


44 


48 


62 


73 


107 


97 


1 


10 


20 


21 


30 


39 


39 


48 


43 


10 


9 


30 


35 


39 


54 


65 


82 


66 


0 


7 


6 


19 


20 


21 


17 


23 


24 


15 


14 


3 1 


46 


63 


59 


75 


90 


108 


7 


9 


25 


42 


35 


60 


59 


5 C 


87 


8 


13 


15 


29 


27 


26 


39 


45 


63 


1 1 


20 


34 


50 


47 


74 


79 


82 


103 


1 


5 


11 


11 


1 0 


15 


23 


30 


21 


5 


31 


29 


35 


61 


69 


93 


108 


91 


2 


7 


9 


14 


16 


24 


29 


29 


33 


3 


7 


21 


29 


24 


30 


26 


35 


55 


7 


12 


21 


24 


3 1 


43 


49 


50 


71 


3 


8 


9 


19 


22 


27 


25 


29 


33 


0 


1 


2 


5 


8 


8 


5 


12 


8 


7 


8 


17 


18 


25 


32 


28 


43 


31 


15 


24 


30 


41 


53 


56 


77 


96 


106 


1 


3 


5 


9 


13 


13 


13 


14 


22 


8 


6 


16 


13 


16 


26 


29 


23 


29 


7 


14 


28 


28 


27 


32 


37 


48 


54 


6 


12 


17 


32 


25 


37 


51 


77 


98 



80 



RANDOM NOS 



FROM THE NEYMAN TYPE A DISTRIBUTION 



Ml VALUE 

12 3 



7 


4 


5 


1 1 


7 


3 


15 


12 


7 


1 


4 


5 


7 


2 


4 


8 


7 


4 


5 


1 8 


7 


7 


18 


11 


7 


1 


7 


8 


7 


3 


6 


7 


7 


16 


26 


27 


7 


2 


8 


6 


7 


17 


25 


41 


7 


4 


16 


10 


7 


2 


3 


13 


7 


18 


21 


45 


7 


8 


11 


14 


7 


4 


11 


22 


7 


15 


17 


26 


7 


14 


15 


23 


7 


5 


c 


21 


7 


4 


18 


24 


7 


5 


13 


29 


7 


4 


6 


19 


7 


7 


9 


10 


7 


6 


6 


14 


7 


8 


18 


33 


7 


2 


6 


8 


7 


4 


11 


27 


7 


E 

-> 


6 


13 


7 


6 


13 


20 


7 


1 C 


24 


1 9 


7 


10 


27 


40 


7 


10 


10 


28 


7 


7 


13 


30 


7 


1 


12 


8 


7 


8 


14 


17 


7 


3 


14 


1 5 


7 


6 


12 


23 


7 


6 


18 


28 


7 


5 


22 


31 


7 


2 


9 


7 


7 


6 


7 


15 


7 


5 


2 


27 


7 


4 


8 


17 


7 


1 1 


15 


27 


7 


10 


24 


24 


7 


6 


19 


18 


7 


4 


12 


1 7 


7 


2 


12 


12 


7 


6 


16 


24 


7 


12 


28 


39 



VALUES 



5 


6 


7 


8 


9 


17 


31 


27 


26 


34 


21 


24 


3 1 


43 


40 


11 


17 


18 


12 


23 


12 


15 


18 


17 


16 


30 


30 


33 


30 


37 


24 


37 


34 


46 


61 


18 


23 


18 


21 


25 


10 


18 


1 7 


30 


23 


53 


68 


85 


84 


100 


14 


20 


lo 


27 


23 


69 


75 


92 


111 


109 


30 


41 


33 


45 


49 


23 


23 


31 


39 


39 


75 


1 12 


89 


92 


121 


30 


32 


28 


48 


66 


35 


31 


33 


48 


54 


50 


55 


76 


74 


91 


42 


55 


64 


76 


74 


32 


30 


53 


45 


53 


42 


46 


55 


62 


64 


35 


48 


70 


65 


41 


21 


24 


26 


23 


36 


36 


36 


43 


45 


52 


23 


37 


43 


45 


57 


45 


48 


57 


71 


87 


8 


16 


24 


18 


31 


40 


51 


54 


73 


50 


20 


20 


31 


35 


39 


39 


32 


46 


52 


70 


54 


65 


6 6 


84 


83 


60 


58 


34 


87 


94 


37 


44 


50 


55 


81 


49 


57 


59 


57 


77 


20 


27 


3 1 


29 


29 


21 


31 


43 


56 


49 


38 


30 


45 


42 


55 


43 


47 


49 


57 


77 


36 


40 


52 


48 


69 


45 


58 


75 


77 


89 


17 


30 


25 


28 


34 


25 


19 


36 


34 


42 


25 


25 


40 


40 


37 


12 


22 


38 


28 


39 


49 


50 


53 


77 


78 


42 


58 


63 


76 


81 


26 


42 


48 


67 


48 


32 


45 


40 


47 


40 


25 


27 


26 


37 


38 


45 


57 


83 


76 


80 


58 


33 


83 


96 


98 



M2 

4 

18 

19 

13 

11 

20 

27 

9 

12 

34 

10 

51 

19 

23 

61 

27 

17 

42 

40 

28 

33 

31 

14 

26 

22 

51 

12 

29 

20 

24 

38 

44 

39 

39 

22 

22 

19 

28 

30 

46 

16 

14 

26 

15 

34 

37 

16 

26 

10 

41 

43 



81 



8 

8 

8 

8 

e 

8 

8 

8 

8 

8 

8 

3 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

3 

8 

a 

8 

8 

3 

a 

8 

8 

8 

3 

3 

6 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 



MCS. FROM THE N EYM AN TY°E A DISTRIBUTION 



I 2 3 



18 


20 


32 


5 


10 


7 


17 


30 


52 


6 


20 


24 


6 


17 


35 


9 


20 


34 


5 


13 


20 


5 


20 


25 


17 


18 


35 


1 0 


8 


39 


7 


21 


28 


16 


28 


47 


6 


18 


21 


5 


11 


25 


2 


17 


26 


7 


12 


25 


8 


21 


41 


10 


12 


24 


2 


4 


7 


8 


14 


9 


4 


17 


16 


12 


18 


26 


8 


20 


23 


5 


21 


35 


13 


25 


46 


Q 


12 


2 1 


7 


8 


16 


8 


21 


30 



6 


8 


8 


12 


18 


32 


7 


18 


26 


3 


10 


4 


5 


16 


25 


11 


29 


45 


7 


25 


26 


7 


20 


24 


7 


9 


13 


1 G 


22 


22 


13 


23 


30 


7 


30 


38 


5 


9 


9 


8 


11 


14 


5 


10 


1 1 


13 


16 


25 


9 


14 


2 2 


10 


18 


29 


15 


23 


39 


7 


22 


41 


9 


18 


27 


7 


13 


21 



8 2 VALUES 



4 



42 

18 

72 

26 

26 

40 
21 
27 

27 
36 
25 
46 
38 

23 
22 

25 

41 
36 

7 

24 

28 
20 
35 
48 
56 

27 
22 
41 
10 

28 

26 
9 

44 
51 
47 
28 
33 
32 

45 
47 
13 
24 

9 

27 



52 

42 

22 

22 



5 


6 


7 


8 


9 


4 6 


57 


65 


83 


99 


19 


30 


3 6 


31 


33 


77 


87 


10 2 


112 


13 6 


33 


42 


46 


53 


82 


41 


54 


76 


76 


80 


53 


54 


65 


74 


92 


34 


28 


34 


53 


48 


44 


48 


67 


75 


89 


55 


71 


77 


fc £ 


98 


54 


52 


68 


69 


7 1 


45 


35 


50 


65 


70 


65 


77 


97 


1 17 


125 


32 


43 


5 1 


57 


59 


37 


51 


62 


70 


86 


40 


46 


46 


67 


81 


3 8 


52 


52 


55 


6 1 


41 


61 


74 


95 


96 


47 


65 


62 


9 1 


107 


1 9 


18 


25 


31 


20 


32 


34 


52 


48 


50 


28 


24 


31 


26 


56 


5 8 


74 


66 


87 


81 


42 


44 


43 


75 


82 


75 


61 


87 


99 


100 


85 


1 12 


10 1 


89 


144 


28 


36 


35 


5 1 


58 


34 


32 


41 


42 


51 


56 


56 


78 


92 


95 


25 


28 


37 


36 


32 


45 


60 


62 


57 


86 


47 


65 


63 


78 


81 


1 1 


26 


38 


34 


30 


48 


55 


65 


81 


102 


63 


67 


89 


1 16 


103 


57 


41 


6 6 


62 


100 


31 


44 


45 


66 


63 


30 


34 


49 


5 c 


67 


44 


54 


58 


66 


72 


3 5 


59 


61 


75 


93 


5 1 


49 


63 


72 


92 


1 1 


19 


15 


17 


31 


28 


35 


48 


53 


53 


25 


37 


29 


40 


35 


50 


48 


67 


70 


73 


25 


52 


50 


52 


57 


58 


60 


82 


74 


103 


67 


93 


73 


82 


130 


6 5 


80 


94 


97 


103 


40 


54 


6 1 


65 


33 


3 C 


23 


41 


44 


5 0 



82 



s 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 



NOS. FROM THE NEYMAN TYPE A DISTRIBUTION 



M2 VALUES 



1 


2 


3 


4 


5 


6 


7 


8 


9 


1 1 


12 


13 


17 


28 


29 


24 


31 


53 


1 1 


21 


24 


33 


4 C 


47 


48 


74 


92 


6 


20 


22 


31 


37 


50 


46 


63 


63 


8 


12 


19 


20 


23 


30 


39 


59 


61 


20 


35 


52 


70 


75 


1C4 


98 


137 


156 


20 


23 


34 


62 


73 


32 


88 


98 


110 


6 


21 


25 


32 


50 


53 


63 


63 


63 


8 


14 


29 


35 


47 


41 


6 0 


67 


67 


A 


9 


14 


22 


24 


34 


27 


37 


43 


12 


17 


36 


45 


61 


83 


74 


89 


87 


9 


20 


30 


57 


67 


51 


75 


1C4 


123 


13 


19 


23 


31 


32 


46 


55 


61 


71 


12 


18 


39 


51 


66 


80 


83 


102 


109 


8 


18 


23 


37 


38 


39 


60 


76 


73 


15 


20 


38 
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